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ABSTRACT 


The ability to use flight data to determine an aircraft model with structural dynamic effects suitable for piloted 
simulation and handling qualities analysis has been developed. This technique was demonstrated using SR-71 
flight test data. For the SR-71 aircraft, the most significant structural response is the longitudinal first-bending 
mode. This mode was modeled as a second-order system, and the other higher order modes were modeled as a time 
delay The distribution of the modal response at various fuselage locations was developed using a uniform beam 
solution, which can be calibrated using flight data. This approach was compared to the mode shape obtained from 
the ground vibration test, and the general form of the uniform beam solution was found to be a good representation 
of the mode shape in the areas of interest. To calibrate the solution, pitch-rate and normal-acceleration instrumenta- 
tion is required for at least two locations. With the resulting structural model incorporated into the simulation, a 
good representation of the flight characteristics was provided for handling qualities analysis and piloted simulation. 

INTRODUCTION 

Performance requirements generally call for a structure to be as light as possible. For large cruise vehicles, 
such as the high-speed civil transport and single-stage-to-orbit vehicles, maneuvering requirements are generally 
low. The combination of low maneuvering requirements with a high vehicle gross weight often results in a relative- 
ly flexible structure where the frequencies of the structural modes approach those of the rigid-body modes. There- 
fore, it is highly desirable to include the structural modes in the analysis of the handling qualities and the piloted 
simulation of the vehicle. A simulator study of a vehicle that included structural dynamic modes concluded that 
structural dynamics can have a significant effect on the handling qualities of the vehicle. However, most piloted 
simulations do not include structural dynamic modes primarily because of the complexity of the structural modes. 
For a full-envelope simulation, the structural mode calculations can be of a similar order of magnitude as those for 
the rigid-body modes, which can often tax the capability of the simulation computers. 

Another difficulty arises when providing simulation support for flight research programs. In this case, the 
structural modes that are available from either a prediction or ground vibration test (GVT) often do not correspond 
to the current configuration being tested. It is generally not feasible from a cost and time standpoint to update the 
structural dynamic models to support handling qualities experiments. The goal of this study was to develop an easy 
method to implement a structural model that can be calibrated with flight data. Because this model is being devel- 
oped to support piloted simulation and handling qualities analysis, it was preferred to have comparatively fewer in- 
strumentation and flight test maneuvers than for the handling qualities experiment itself. This report documents the 
development and calibration of the structural models using actual flight test data. 

In the SR-71 program, the flight data have shown a noticeable structural response. The structural modes have 
generally not been modeled in the simulations or past handling qualities analyses. Many handling qualities cntena 
today depend on frequency response characteristics that generally extend to frequencies that include some structur- 
al dynamic mode effects. As a result, it was necessary to expand the current SR-71 models to include structural dy- 
namic effects. The first fuselage bending dynamics were easy to distinguish. Although the pilots were very 
conscious of the structural dynamic motion, it was not a handling qualities problem because the motion was far 
enough removed from the rigid-body frequencies. The mode was well enough defined and of sufficient magnitude 
that it could be used to develop the modeling method for incorporating the structural dynamic effects. The SR-71 
aircraft was used to validate the simulation and analysis capability. This modeling technique provides the capabili- 
ty to use flight data to establish structural models suitable for simulation and analytical studies that will be applica- 
ble to other large, flexible vehicles of the future. 



NOMENCLATURE 


A normal acceleration, g 

n 

Aj mode shape bias 

A 2 mode shape angular scaling factor 

c.g. center of gravity, percent 

C r control surface input effectiveness coefficient, per deg 

6 

dz/dx slope of normalized structural mode vertical deflection 
F g control surface input effectiveness, in/deg 

FS fuselage station, in. 

FSq reference fuselage station, in. 

g acceleration caused by gravity, ft/sec 2 

GVT ground vibration test 

Kj displacement constant 

K 2 slope constant 

L characteristic length, ft 

L a rigid-body lift caused by angle of attack, 1/sec 

Fg rigid-body lift caused by elevator, 1/sec 

M Mach number 

M rigid-body pitching acceleration caused by pitch rate, 1/sec 

M a rigid-body pitching acceleration caused by angle of attack, 1/sec 

Mg rigid-body pitching acceleration caused by elevator, 1/sec 

q pitch angular rate, deg/sec 

q dynamic pressure, lb/ft 2 

q pitch angular acceleration, deg/sec 2 

s Laplace operator 

S wing area, ft 2 

SAS stability augmentation system 

V velocity, ft/sec 

W aircraft gross weight, lb 

x horizontal axis coordinate, ft 

z vertical displacement, in. 

z vertical velocity, in/sec 

8 elevator deflection, deg 

5 p pilot stick input, in. 

Ax distance from center of gravity to sensor, ft 

A8 incremental elevator deflection from trim, deg 
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^ damping ratio 

tj structural modal displacement, in. 

f) structural modal velocity, in/sec 

tj structural modal acceleration, in/sec 2 

0 pitch angle, deg 

x normal-acceleration time delay, sec 

a n 

x 6 elevator time delay, sec 

to natural frequency, rad/sec 

Subscripts 

aft aft instrumentation location 

d delayed 

fwd forward instrumentation location 

r rigid value 

s structural value 

0 bias term 

STRUCTURAL DYNAMIC MODEL 

The most significant structural dynamic response of the SR-71 aircraft is the longitudinal first-bending mode. 
Figure 1 shows the flight-measured pitch rate and normal acceleration caused by a sharp, longitudinal pitch stick 
input as well as the rigid-body motion obtained from the simulation response to the same input. The flight response 



Figure 1 . Comparison of flight response at forward-sensor location with rigid-body simulation response. 
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has a considerable additional higher frequency component because of the structural response. A structural response 
of this type can be described by the product of the dynamic response of the mode and the spatial distribution of the 
structural mode (i.e., the mode shape): 

Z = T| Kj(x) 

The dynamic portion of the response, T|, is the response at one point along the fuselage, usually of maximum 
amplitude. The mode shape, Kj (x), is the distribution of the response for other locations along the fuselage. The 
exact distribution of the deflection can be obtained from a structural analysis using the actual mass and structural 
characteristics from a GVT or from flight data. In general, the distribution would be a function of loading and flight 
condition. In the following section, an approximate distribution will be developed based on the solution for a uni- 
form beam. It will be shown that this distribution provides a reasonable representation of the structural mode for 
the purpose of performing handling qualities analysis. Flight data can be used to calibrate the approximate solution 
with a minimum of instrumentation to provide the best match of the actual vehicle as tested. 

Mode Shape Model 

The displacement constant, K,, is the distribution of the response for locations along the fuselage. For a uni- 
formly loaded beam, the displacement constant is a sinusoidal function: 

Kj = sin(p x/L + 0) + Aj 

The nodes of the first-bending mode are located at x = 0.224 L and at x = 0.776L, where L is the beam 
length. 2 The other constants are p = 1.57c, <|) = 0.7571, and A, = 0.267. The primary variable is the nondimen- 
sional beam length x/L . Figure 2 illustrates the resulting mode shape. 
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The uniform beam solution has three variables that can be used to fit a particular mode shape: L, p, and Aj. 
The form of the solution can be rearranged slightly to provide more meaningful parameters. Because the structural 
coordinate system of the uniform beam solution is not necessarily coincident with the aircraft coordinate system, a 
bias between the two coordinate systems is required so that x = (FS - FS Q )/ 12, where FS is the fuselage station 
of the aircraft reference system in inches. For easier interpretation of the variables, it is useful to change the phase 
parameter so that FS 0 corresponds to the center of the beam where the slope is zero. The resulting expression for 
the deflection constant is as follows: 


K, = -cos 


1.5tc 


FS-FS 0 > 
12L 


+ Ai 


To test the validity of the assumption that a uniformly loaded beam solution can reasonably represent an actual 
aircraft mode shape, this form of solution was used to fit the mode shape obtained from a GVT. The GVT data ob- 
tained from the YF-12 aircraft, 3 which is geometrically similar to the current SR-71 test aircraft, provide a defini- 
tion of the shape of the first-bending mode. A curve fit of the data was obtained using the sinusoidal form of the 
uniform beam solution. Figure 3 shows the GVT data and the curve fit. The fit of the GVT data results in an effec- 
tive length of 139 ft (compared to the actual fuselage length of 107 ft) and an inflection point location of FS 0 798. 
The bias, Aj, was 0.72 rather than the uniform beam solution of 0.267. The uniform beam form of solution pro- 
vides a good fit of the actual aircraft deflection curve from the GVT data, except near the aft end of the airplane. 
For handling qualities analysis, the range of interest is from the nose to the center of gravity ( c.g .) location 
(FS 900), which encompasses the pilot location, control system, and instrumentation sensors. 


Normalized 

deflection 



5 



Dynamic Response Model 


The dynamic portion of the response, T|, is modeled as a second-order system with the elevator as the forcing 
function: 


ti + 2£(ofi + (OT| = G) F 5 A8 


The term F 6 is the effectiveness of the elevator which, in conjunction with the normalized deflection K,, 
scales the structural deflection response z s in inches when the elevator input is in degrees. The pitch angular de- 
flection of the structure is the slope of the deflection curve with respect to the longitudinal distance and is given by 
the following: 

dz/dx = — K 2 t] 

The slope dz/dx is approximately the angle of the fuselage relative to the rigid-body x-axis.* The slope con- 
stant, K 2 , is given by the following: 


-A, 


k 2 = 


■sin 


1.5rc 


FS-FS 0 > 

\2L 


The scale factor, A 2 , provides the proper scaling for the pitch angle equations. The structural dynamic dis- 
placements and angles are perturbations from the rigid-body values. The structural dynamic components of the 
various response parameters are as follows: 


Displacements: 

z s = Kj T1 

(in.) 

CD 

II 

* 

to 

(deg) 

Rates: 

z s = K , n 

(in/sec) 

«s = K 2 *1 

(deg/sec) 

Accelerations: 

K l .. 
" 12^ 

(g) 

= k 2 Ti 

(deg/sec 2 ) 


For the nonlinear simulation, the elevator input to the structural mode, A8, is the incremental elevator deflec- 
tion from trim. The effects of static structural deflections are included in the basic aerodynamic model, which con- 
tains the static aeroelastic corrections. The complete equations for the two types of measured parameters that were 
available on the test aircraft, normal acceleration and pitch rate, are as follows: 


A = 7 ^— 
"» 12 g 


-cos 


1.571 


(FS-FS o y 


12 L 


+ A 1 — 

. Ls + 


2 2 
S (0 


2a>£s + a> J 


F S A8 


r -a 


= 


sin 


1.57t 


(F5-F5 0 ) > 
12L 


s (0 


2 2 
_s + 2 cd£s + co . 


F 6 A5 


The structural dynamics are then added to the rigid-body motion to obtain the total vehicle response. To ac- 
count for the effects of higher order structural modes that have not been modeled, a time delay was introduced be- 
tween the elevator deflection and the input to both the rigid-body aerodynamics and the structural modes. This 


♦Because of the sign convention of x (+ x toward the rear) and z(+ z up), a positive slope corresponds to a negative pitch angle. 
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procedure is similar to that used in handling qualities analysis to determine lower order equivalent systems where 
the higher order control-system dynamics are represented as a time delay. The plausibility of such a delay can be 
visualized by considering an aircraft that has a wing with torsional flexibility. The initial result of a surface deflec- 
tion will be a wing twisting. After the wing has twisted, aerodynamic forces will be generated, resulting in rigid- 
body pitch acceleration and longitudinal first-bending mode acceleration. This phenomenon was represented as a 
delay between the surface position as produced by the actuator and the surface position that is used to drive both 
the rigid and flexible responses. An additional delay was allowed for the acceleration response as shown in the fol- 
lowing transfer functions: 


8 d (s) -T g s 

8(s) ~ ‘ 

q(s) _ 

'q r (s) qJsY 

5 d (s) 

8 d (s) + S d (s) 

A n (s) 

" A nc« r (s) Ax 

s d < s > 

8 d (s) + g 


tf r ( s ) 

s d (s) 


+ 



S 


For the current SR-71 flight test program, there was no instrumentation for the control-surface positions. The 
control-surface responses were calculated using the simulation models for the actuators and the pilot control input, 
which was recorded. As a result, it was not possible to discriminate between delays caused by structural effects and 
delays caused by errors in the actuator models. 


CALIBRATION TECHNIQUES USING TEST DATA 


The unknown variables required to define the structural model are <o s , £ s , FS# U A,, A 2 , F g , X , and T g . 
Flight data were available for normal acceleration and pitch rate at two fuselage stations: one at FS 235, near the 
nose and the other at FS 683, near the c.g. The best method to solve the unknown parameters is a parameter- 
estimation technique such as the pEst software, 5 which readily identify the unknown structural characteristics and 
the generally unknown rigid-body aircraft characteristics. The time delay can be estimated by shifting the input 
time history various amounts 6 or solved directly by using the nonlinear estimation capability of the pEst program. 
For the SR-71 aircraft, the structural mode identification can also be accomplished with a simplified analysis be- 
cause the structural mode frequency is sufficiently removed from the short-period frequency and the stability aug- 
mentation system is not significantly coupled with the structural mode because of the sensor location. In general, 
these conditions will not exist. The motivating factor for including the structural modes in the handling qualities 
analysis is that the structural frequency is near enough to the short-period frequency to cause a problem for pilot 
control. This also makes a simplified analysis difficult and generally can only be analyzed by using parameter- 
estimation techniques. The following section will develop the method used for the simplified analysis to provide 
some insight into the results. 


Simplified Analysis Technique 

The four mode shape parameters— FS 0 ,L, A,, and A 2 — can be determined from the ratios of the magnitudes 
of the two pitch-rate and two normal-acceleration responses. These ratios are A n / q at the forward and aft loca- 
tions, A /A and <? fwd /<? aft • When these ratios are obtained from the free-oscillation part of the response fol- 
lowing a sharp excitation input, the response is primarily caused by the structure. From the previous expressions 
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for A n and q s , the ratios can be evaluated at the structural frequency co s . The phase information can be used to 
determine the location of the motion along the x-axis. When the positive A n peak leads the positive q s peak by 90°, 


the ratio 


is considered positive. The ratios 


fwd 


*aft 


and 


’fwd 


’•ft 


are considered positive when both responses 


are on the same side of the node. Figure 4 provides an example of the free-oscillation portion of a time history. In 
this example, the amplitude ratios can be read directly from the time history. However, in many cases, it was more 
difficult to directly measure the amplitudes because of the noise on the measurements. As a result, a curve fit of 
each parameter was made and the amplitudes were determined from the fitted curves. 


*n> 

9 



q. 

deg/sec 


Figure 4. Free-oscillation portion of response to sharp pitch control input, M 0.8. 


For the example time history, the measured value of the structural mode frequency was 14.6 rad/sec. The am- 
plitude ratios from the measurements at the forward location (FS 234.5) and the aft location (FS 683) resulted in 
the following equations: 

• Ratio of acceleration to pitch rate at the two fuselage stations 


cos 


1.5ti 


(234.5 - FS q )' 

12Z 


-A, 


fwd 


12gA 2 


sin 


cos 


1.571 


c _(234.5 -FS 0 ) 
1-5 12L 

(683 -FS 0 )) 

12L 


14.6 = 0.058 g/ deg/sec 


-A, 


aft 


12gA, 


sm 


1.57t 


(683 - FS 0 f 

HL 


14.6 = -0.034 g/ deg/sec 
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• Ratio of accelerations at forward and aft fuselage stations 


fwd 




cos 


1.5tt 


(234.5 -FS 0 ) 

\2L 


-A, 


= -5.72 


COS 1.571 


(683 -FS 0 ) 
12 L 


-A, 


• Ratio of pitch rates at the forward and aft fuselage stations 


5 fwd 


51 aft 


sin 1.5 n 


(234.5 -FS 0 )\ 


12L 


) _ 


3.40 


sin 


1.5 Jt 


(683 - FS 0 ) 


12 L 


These equations were then solved for the four unknowns L, FS 0 , A p and A 2 . The solution required an itera- 
tive procedure, and a spreadsheet with an equation solver was used for this process. Good starting values for the so- 
lution are L = fuselage length and FS 0 = L/Z Figure 5 shows the solution. It can be seen that the assumption of 
the sinusoidal form for the deflection distribution provides a rather nonlinear solution obtained from a minimum 
number of measurements. The next step was to determine the structural gain F g and the time delays. This was ac- 
complished by putting the mode shape data into the simulation and then adjusting the gain and delays to obtain the 
best fit to the flight data. The delays were adjusted only to the nearest simulation frame (0.020 sec). The delay in 
the elevator control input, determined from the pitch-rate response, was found to be 3 frames or 0.060 sec. A 
4-frame delay (0.080 sec) had to be added to the normal-acceleration responses. This is likely caused by delays in 
the accelerometer instrumentation itself rather than any structural phenomena. 



Figure 5. Pitch-rate and normal-acceleration amplitude ratios and calculated fit. 
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Parameter-Estimation Techniques 


Two parameter-estimation techniques were used: the pEst software 5 and the spreadsheet identification method, 
which uses the equation solving capability of a spreadsheet on a personal computer. These techniques use the lon- 
gitudinal short-period rigid-body equations of motion and structural equations of motion. The input to the system 
of equations was the elevator deflection. Because the elevator deflection was not measured in flight, a calculated 
value was obtained by using the measured stick position and the simulation to obtain the elevator position. The el- 
evator position was the average of the inboard and outboard elevator positions. The simulation actuator models 
were second order, with the inboard elevator driven by the stick command and the outboard elevator driven by the 
inboard elevator position. The elevator input to the aerodynamic and structural models was delayed by x g so that 
the input at time t was 8(r - x g ). The following equations were used to model the dynamic responses: 

q r (t ) = M q q r (t) + M a a r (t) + M g 8(r - x g ) + M 0 
a r (0 = q r i t)-L a a r (r)-L g 8(/-x g ) + L 0 
ri(r) = -2 £ to ii(r) -co 2 ti(r) + © 2 F g 6(/-x g ) + F 0 
The output equations were as follows: 


v( ( 


Ax . ( 

^ i 

f 

= ;K 

-\j-4-\jj 


a J Mg 

-cos 

V 


1.571 


(FS-FS^ 

ML 


+ A, 


fi[r-x„ ) + A„ 
\ V "o 


<?(0 = 9 r (0 + 


-a 2 . ' 

— r-sm 


1.57X 


( FS-FS 0 ) 

ML 


Yi 


11(0 + <?o 


where FS and Ax were evaluated for the forward and aft instrument locations. 


The cost function was the square of the difference between the calculated outputs and the flight-measured val- 
ues of the four measurements. The unknown variables used to minimize the cost function were the structural coef- 
ficients, short-period rigid-body coefficients, and various bias terms: a), £, F g , L, FS 0 , Aj, A 2 , x g , x fl , M q , M a , 

M h' L a' L & V ’ L 0’ M 0> F 0- A n 0 ’ A n 0 ’ ' 

fwd aft 

For the spreadsheet solution, the time delay was determined by evaluating the cost function for various time 
delays, which were incremented as multiple sample intervals. 6 This was done for only one flight condition, and the 
results were used for the other conditions. For the pEst technique, the time delays were implemented as interpolat- 
ed values of tables (i.e., 8 as a function of time) and the time delay was estimated directly for each time history. 
With both techniques, the parameters of the nonlinear mode shape could be solved for directly. Other parameter- 
estimation techniques that require linear equations could be used by estimating Kj and K 2 for the forward and aft 
location. The mode shape parameters L, FS 0 , Aj, A 2 could then be calculated from K, and K 2 . 


Figure 6 provides a comparison of the mode shapes obtained from the simplified analysis results and the 
parameter-estimation technique results. All of the results show good agreement in the area of interest from the nose 
to the c.g. Some disagreement is seen at the aft end of the aircraft; however, if this area were of importance, 
additional instrumentation could be added in this region. Figure 7 shows a typical example of the fit of the flight 
data using the spreadsheet parameter-estimation technique for a flight condition at M 3.0. A reasonably good fit 
was generally obtained for all the flight conditions. 
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Figure 6. Comparison of mode shapes resulting from three analysis techniques, M 0.8. 


SUMMARY OF PARAMETER-ESTIMATION RESULTS 

The data used to analyze the structural characteristic data were collected on five flights of the SR-71 airplane. 
The maneuvers consisted of sharp raps of the longitudinal control stick, which primarily excited the structural 
mode In one case, a frequency sweep maneuver was performed. A summary of the flight conditions and maneu- 
vers is contained in table 1. The maneuvers at M 3.0 were performed with the stability augmentation on; the other 
maneuvers were performed with the stability augmentation off. The following discussion will present the results 
obtained from the five flights for the simplified analysis, spreadsheet parameter identification analysis, and pEst pa- 
rameter identification analysis. 


Table 1. Flight conditions for the test maneuvers. 


Right 

Maneuver 

Mach 

number 

Altitude, ft 

Dynamic 

pressure, 

lb/ft 2 

Aircraft 
weight, lb 

15 

Stick rap 1 

0.785 

20,200 

440 

72,900 

15 

Stick rap 2 

0.800 

19,800 

440 

72,400 

18 

Stick rap 1 

0.825 

18,000 

490 

120,000 

18 

Stick rap 2 

0.825 

18,400 

490 

119,900 

18 

Stick rap 3 

0.825 

18,700 

490 

119,800 

19 

Stick rap 1 

0.850 

22,800 

400 

115,900 

19 

Stick rap 2 

0.850 

23,200 

445 

115,800 

20 

Stick rap 1 

0.800 

20,600 

426 

11,750 

20 

Stick rap 2 

0.820 

21,000 

440 

117,100 

20 

Stick rap 3 

0.830 

21,200 

447 

116,700 

21 

Stick rap 1 

3.000 

75,300 

460 

80,300 

21 

Frequency sweep 

3.000 

75,300 

460 

80,000 
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Figure 7. Time history of match between flight data and calculated time history using spreadsheet identification 
technique, M 3.0. 
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Mode Shape 

Figure 8 shows the results of all the points in terms of the mode shape parameters L, FS Q , A v and A 2 . The 
simplified analysis and the two parameter identification technique results showed good agreement, although the pa- 
rameter identification results were usually more consistent. Good agreement was found between the spreadsheet 



Weight, lb 950294 

Figure 8. Mode shape parameters L, A,, FS Q , and A 2 as function of aircraft gross weight, (open symbol, M 0.8; 
solid symbol, M 3.0). 
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parameter-estimation technique and the pEst technique. The faired line shown in figure 8 is a least square fit of the 
two sets of parameter identification data. The parameters FS 0 and A, were the most significant in determining the 
mode shape. Figure 9 shows the effect of the parameter variations with weight on the mode shape for two weight 
conditions. A fairly large shift in the forward node toward the aft is seen as the gross weight is increased. 


Normalized 

daflaction 



Figure 9. Comparison of mode shapes for heavy-weight and lightweight aircraft using faired data from figure 8. 


Control Effectiveness 

Figure 10 shows the control effectiveness parameter. The control effectiveness was modeled as Fg = Cf & . It 

would be expected that the control effectiveness for the structural mode would be similar to the control effective- 
ness for the rigid-body motion. A curve proportional to pitching moment caused by elevator is shown in figure 10. 
The flight data are in general agreement with this trend although data as a function of Mach number are very 



Figure 10. Structural control effectiveness parameter. 
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limited. A simple approximation to this function of Mach number is a function inversely proportional to Mach 
number with an upper limit, as shown in figure 10. 


Frequency and Damping 


Figure 11 shows the structural frequency and damping. The frequency was a function of vehicle gross weigh 
but was not significantly affected by flight condition. The damping ratio was independent of airplane gross _ weight 
but was noticeably reduced at the M 3.0 flight condition. Figure 12 shows the decrease in damping with Mach 
number A function similar to the one used for the control effectiveness is also shown as an interpolation over the 



Figure 11. Structural frequency and damping characteristics as function of airplane weight (open symbol, M 0.8; 
solid symbol, M 3.0). 



Figure 12. Structural damping characteristics as function of Mach number. 
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Mach number range; however, there are only marginal data to support this trend. No effect on damping was seen 
because of dynamic pressure over the limited range tested. 

Time Delay 

The time delay was estimated with the spreadsheet parameter-estimation technique at one flight condition. 
Several fits were obtained using a range of time delays for both the elevator and the normal acceleration. Table 2 
lists the cost function results from the matches. An interpolation of these data produced an estimate for the elevator 
delay of T g = 0.064 sec and t a = 0.077 sec. It can be seen that the simplified analysis estimate of a 3-frame de- 
lay (0.060 sec) for the elevator and a 4-frame delay (0.080 sec) for the normal acceleration agrees well with the 
spreadsheet parameter-estimation results at the same flight condition. An easier estimation of the time delay was 
obtained from the pEst program, which had the capability to directly estimate the time delays for each maneuver. 
Figure 13 shows the three sets of results. Considerably more scatter is seen when all of the pEst points are includ- 
ed. The pEst results indicated an average elevator time delay of t g = 0.036 sec and an average normal- 
acceleration time delay of x a =0.100 sec. 


Table 2. Cost function from the spreadsheet parameter estimation as a function 
of pitch-rate and normal-acceleration time delays. 


.see 

n 




0.06 

0.08 

0.10 


0.04 

18.5771 

16.2385 

20.1536 

x 5 , sec 

0.06 

12.2156 

10.696 

13.9351 


0.08 

14.7924 

13.6397 

16.7074 
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Average 
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(a) Elevator. 

Figure 13. Elevator and normal-acceleration time delays. 
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(b) Normal acceleration. 

Figure 13. Concluded. 
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SIMULATION RESULTS 


Simulation Model 

A simulation model of the longitudinal first-bending mode was created from the parameter estimation results. 
The structural mode equations were added to the rigid-body states. A low-pass filter was used to approximate the 
steady-state trim elevator deflection. The incremental elevator deflection, which is the input to the structural mode 
model, is the difference between the current and the filtered elevator positions. The transfer functions for the struc- 
tural motion and the combination with the rigid-body motion are shown in the following equations as a function of 
the fuselage station: 




V s > 


1 


A8(s) 12g 


-cos 


1.5n- 


FS - FS. 


\2L 


<?s( s ) 

A8(s) 


-An 


: sm 


( FS-FS, 
1.5n- 


12L 


HI; ris: 

JLs + 2<d£s + (0 . 


(0 


*,(»> 

5(s) 


C *r 


(S) 


Ax <? r (s) . \ ( * r 


5(S) 


+ — 
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<7 r (s) _ 9 s (s)]_-t 8 s 
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The control effectiveness was modeled as a function of Mach number and dynamic pressure as follows: 



= 

0.06 IM -0.01 

c h 

< 

0.065 


= 

&C F 
W F 6 

The structural mode shape parameters were modeled as a function of airplane gross weight as follows: 

L = 

208 

- 0.000105 W 

FS 0 = 

752 + 0.00063 W 

A, = 

0.803 + 0.00000049 W 

a 2 = 

113 

- 0.000214 W 


As shown below, the structural frequency was modeled as a function of airplane gross weight, and the structural 
damping ratio was varied as a function of Mach number: 

co s = 23.6 - 0.000079 V/ 

C s = 0.075 IM- 0.007 

C s < 0.087 

The time delays were modeled as pure time delays with the following values: 

Tg = 0.036 sec 

= 0.100 sec 

n 

Comparison of Flight/Simulation Frequency Response 

The frequency response characteristics are of particular interest for handling qualities analysis. A frequency 
sweep maneuver was performed in flight at M 3.0, and this was used to compare the results from the simulation 
model that was developed with the flight data. The stick input from the flight maneuver was used to produce a sim- 
ulation time history of the frequency sweep maneuver. Fast-Fourier transforms were used on both the flight and 
simulated data to find the pitch rate and normal acceleration to pilot input transfer functions. Figure 14 shows the 
results for the forward instrumentation location, and figure 15 shows the results for the aft instrumentation loca- 
tion. The simulation model shows reasonably good agreement with the flight data at the M 3.0 flight condition. The 
simulation model also allows the extrapolation to other fuselage locations. Figure 16 shows the frequency response 
characteristics for the pilot location, which are of interest for handling qualities analysis. Figure 16 also shows the 
original rigid-body responses from the simulation. Significant improvements in the model fidelity are shown in the 
higher frequency region, especially for the pitch-rate response. The results indicate that the desired objective, to 
provide a model suitable for handling qualities analysis which includes the structural dynamic effects, has been 
achieved. 
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(b) Pitch rate for pilot stick input. 

Figure 14. Comparison of normal-acceleration and pitch-rate transfer functions from flight and simulation for 
forward-sensor location at M 3.0. 
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CONCLUDING REMARKS 


The ability to use flight data to determine an aircraft model with structural dynamic effects suitable for piloted 
simulation and the analysis of handling qualities and flight control system has been developed. This technique was 
demonstrated using SR-71 flight test data. For the SR-71 aircraft, the most significant structural response is the lon- 
gitudinal first-bending mode. This mode was modeled as a second-order system. The effect of other higher order 
modes was modeled as a time delay to both the rigid-body and structural motion. The distribution of the modal re- 
sponse at various fuselage locations was developed using a uniform beam solution, which can be calibrated using 
flight data. This approach was compared with the mode shape obtained from GVT, and the general form of the uni- 
form beam solution was found to be a good representation of the mode shape in the areas of interest. A simplified 
analysis calibration technique was developed and is applicable when structural dynamic motion is not significantly 
altered by the control system and is separated significantly from the rigid-body motion. For the more difficult case, 
where the control-system or short-period interactions prevent a simple analysis, a parameter-estimation technique 
was demonstrated. This analysis can be accomplished by using either standard parameter-estimation programs or a 
personal computer spreadsheet analysis with equation solving capability. To calibrate the uniform beam solution, 
pitch-rate and normal-acceleration instrumentation is required for at least two locations. This technique provided 
structural mode shape information that is comparable to a GVT with a minimum of instrumentation, flight time, 
and analysis. With the resulting structural model incorporated into the simulation, a good representation of the 
flight characteristics was provided for handling qualities and control-system analyses. 
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